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Abstract 

The present state of the nn transition problem is briefly outhned. The models based on 
the diagram technique for direct reactions, potential description of n-medium interaction 
and fleld-theoretical approach are considered. It is shown that for the nn transition in 
medium fleld-theoretical approach should be used. The lower limit on the free-space nn 
oscillation time Tmin is found to be: 10^^ yr > Tmm > 1-2 • 10^ s. 
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1 Introduction 



At present several models of nn transitions in medium are treated. The part of them gives 
radically different results. This is because the process under study is extremely sensitive to the 
details of the model and so we focus on the physics of the problem. 

We briefly outline the present state of the nn transition problem. In the standard calcula- 
tions of ab oscillations in the medium [1-3] the interaction of particles a and b with the matter 
is described by the potentials Uafi (potential model). Imf/5 is responsible for loss of 6-particle 
intensity. In particular, this model is used for the nn transitions in a medium [4-11] followed 
by annihilation: 

n^n^ M, (1) 

here M are the annihilation mesons. 

In [9,10] it was shown that one-particle (potential) model mentioned above does not de- 
scribe the process (1) and thus total neutron- antineutron transition probability: the process 
(1) probability is 1^ ~ F (see Eq. (14)), whereas the potential model gives W ~ l/F [4-11] (F 
is the annihilation width of fi in the medium). In the potential model the effect of final state 
absorption (annihilation) acts in the opposite (wrong) direction, which tends to the additional 
suppression of the nn transition. Since the annihilation is the main effect which defines the 
speed of process (1), the potential model should be rejected. This is because the unitarity 
condition is used for the essentially non-unitary S*- matrix [9,10]. The interaction Hamiltonian 
contains the antineutron optical potential Un and IrnUn plays a crucial role. The ^-matrix 
should be unitary. 

More formally, the expression for the total process width is obtained by means of optical 
theorem. The basic equation J2f^i \ Tfi p~ 2ImTii, S = 1 + iT follows from the unitarity 
condition {SS'^) fi — Sfi. However in the potential model the S-matrix is essenim//y non-unitary 
{SS^)fi = Sfi + afi, afi ^ 0, resulting in X^/^i I Tfi p?« 2ImTii + an ^ llmTn because llmTn 
is extremely small: ^ImTn < 10^^^ [9,10]. The above-given basic equation is inapplicable in 
this case. The potential model describes only the channel with n in the final state [10] when 
unitarity condition is not used. 

For the oscillations in the external field [12,13] the Hamiltonian is hermitian and so there 
is no similar problem. The above-given remark holds only for the processes (1) and total 
neutron- antineutron transition probability calculated by means of non-hermitian Hamiltonian 
(potential model). 

The potential model was developed in 1980-1981. In more recent papers the verious details 
of the model have been refined, in particular the parameters of optical potential. We don't 
dwell on these papers since the heart of the problem is in the non-Hermiticity of the optical 
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potential. For similar reason for the model based on the diagram technique (see below) we 
consider papers [14,15] only. 

In [14] we have proposed the model based on the diagram technique for direct reactions (see 
Fig. 1). This model (later on referred to as the model 1) does not contain the non-hermitian 
operators. For deuteron this calculation was repeated in [15]. However, in [16] it was noted 
that this model is unsuitable for the problem under study. In Sect 2 this problem is considered 
in detail. The model shown in Fig. 2 (model 2) has been proposed [11,16]. 

In the model 2 the antineutron propagator can be bare or dressed. The corresponding 
results differ radically for the following reason. A distinguishing feature of the problem under 
study is the zero momentum transfer in the nn transition vertex. Because of this the (S'-matrix 
amplitudes corresponding to the processes shown in Figs. 2a and 2b with bare propagator 
contain infrared divergence (see Eq. (9)). The small value of the antineutron self-energy E 
removes the singularity from the process amplitude which changes the result radically. This 
circumstance demonstrates the result sensitivity to the details of the model. As a consequence 
the lower limit on the free-space nn oscillation time Tmin given in this paper is in the wide range: 
10^^ yr > Tjnin > 1.2 • 10^ s. In addition, certain of the authors consider that the model 1 (and 
even potential model) is reasonable. From the preceding, it is seen that the problem invites 
further investigation. 

In this paper we compare the model based on diagram technique for direct reactions (section 
2) with the the model based on the field-theoretical approach (section 3) and show that the 
former model is unsuitable for the problem under study. The model 2 and in the first place 
the value of antineutron self-energy invites further investigation. By means of field-theoretical 
approach the range of permissible values of Tmin is obtained (section 3) . Section 4 contains the 
conclusion. 

2 Model 1 

Consider now the model 1 [14]. The Hamiltonian of nn transition is [4] 

Hnn = e^n^n + H.C. (2) 

Here e is a small parameter with e = l/rnn, where Tnn is the free-space nn oscillation time. 

Denote: A = (N.Z), B = {N — 1,Z) arc the initial and intermediate nuclei, M is the 
amplitude of virtual decay A ^ n + (A — 1), M^"^ is the amplitude of nB annihilation in (n) 
mesons. En is the pole neutron binding energy; m, m^, niB are the masses of the nucleon and 
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Figure 1: Model 1 for the nn transition in the nuclei followed by annihilation. 

nuclei A and B, respectively. The amplitude is given by 

(„) ^ JmhriB f M(g)MW(m^) 

For deutron the process probability Wi{t) was found to be 

Wi{t) = ^r^pt, (4) 

where Tnp is the np annihilation width. For carbon and oxygen Wi ~ TnBt/E^ as well (TnB 
is the width of nB annihilation). The lower limit on the free-space nn oscillation time which 
follows from stability of oxygen is [14] 

T^,(O) = l.l-10«s. (5) 

The limit obtained from stability of Fe is Tji,in(Fe) = (8 ^ 11) • 10^ s [15]. 

We list the main drawbacks of the model which are essential for the problem under study: 

1) The model does not reproduce the nn transitions in the medium and vacuum. If the 
neutron binding energy goes to zero, Eq. (4) diverges (see also Eqs. (15) and (17) of Ref. [15]). 

2) Contrary to the model 2 (see next section), the amphtude (3) cannot be obtained from 
the Hamiltonian because in the interaction Hamiltonian there is no term which induces the 
virtual decay (iV, Z) — )■ n + (iV — 1, Z). The neutron of the nucleus is in the bound state and 
so it should be described by the wave function and not the propagator. 

3) The model does not contain the infrared singularity for any process including the nn 
transition, whereas it exists for the processes in the medium and vacuum (see [11,16] and Sect. 
3 of this paper). This brings up the question: Why? The answer is that for the propagator the 
infrared singularity cannot be in principle since the particle is virtual: ^ m? + p^. Due to 
this the model is infrared- free. 
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4) Since the model is formulated in the momentum representation, it docs not describe the 
coordinate-dependence, in particular the loss of particle intensity due to absorption. Also there 
is a no the dependence on nuclear density! The model is crude and has very restricted range 
of applicability. 

We consider the points 2) and 3). The nn transition takes place in the propagator. As 
the result the model is infrared-free. For the processes with zero (or very small) momentum 
transfer this fact is crucial since it changes the functional structure of the amplitude. 

On the other hand, the neutron propagator arises owing to the vertex of virtual decay 
A ^ n + {A — 1). However, as pointed out above, in the interaction Hamiltonian there is no 
term which induces the virtual decay A^ n + (A — l). This vertex is the artificial element of 
the model. It was introduced in order for the neutron (pole particle) state to be separated. 

We assert that for the problem under study this scheme is incorrect. The neutron state is 
described wrongly. The diagram technique for direct reactions has been developed and adapted 
to the direct type reactions. The term "diagram technique for direct reactions" emphasizes this 
circumstance. The processes with non-zero momentum transfer are considerably less sensitive 
to the description of pole particle state. The approach is very handy, useful and simple since 
it is formulated in the momentum representation. This approach was applied by us for the 
calculation of knock-out reactions and p- nuclear annihilation [17]. The price of simplicity is 
that its applicability range is restricted. At the same time, as is seen from p. 3), the process 
under study is extremely sensitive to the description of neutron state. The same is true for the 
value of antineutron self energy (see next section). Besides, the problem is unstable [18]. 

3 Model 2 

Model 2 describes both nn transition in the medium (Fig. 2a) and nuclei (Fig. 2b). It 
corresponds to the standard formulation of the problem: the \in >-state is the eigenfunctions 
of unperturbed Hamiltonian. In the case of the process shown in Fig. 2a, this is the neutron 
plane wave: 

np{x) — Q,"^^'^ exp{—ipx), (6) 

p = (e„,p„), e„ = p^y2m + where f/„ is the neutron potential. In the case of Fig. 2b, 
this is the wave function of bound state [19]. For the nucleus in the initial state we take the 
one-particle shell model. Since the neutron is described by the bound state wave function and 
not the plane wave, the antineutron Green function differs from that in Fig. 2a. To emphasize 
this fact, for the virtual antineutron in Figs. 2b and 2a the different terms are used. 
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Figure 2: Model 2 for the nn transition in the medium (a) and nuclei (b) followed by annihi- 
lation. 



The interaction Hamiltonian is given by 

Hi — 'Hnn + ^) (7) 

where T-L is the hermitian Hamiltonian of n-medium interaction in the case of Fig. 2a and 
Hamiltonian of n-nuclear interaction in the case of Fig. 2b. In so doing the antineutron 
propagator can be bare or dressed. 



3.1 Model 2 with bare propagator 

Consider now the diagram 2a. We show the result sensitivity to the details of the model. If we 
use the general definition of amplitude of antineutron annihilation in the medium Mann which 
is given by 

</| \Onp>= N{2^Y6\pf-pi)Mann, 

= rexp(-i j dxU{x)) - 1 (8) 

(I Oup > is the state of the medium containing the n with the 4- momentum p, N includes the 
normalization factors) then the amplitude corresponding to Fig. 2a diverges: 

Mb = eGbMann: 

Gh = ~ - (9) 

since Pn — Pn, — ^n- This is infrared singularity conditioned by zero momentum transfer 
in the nn transition vertex. Since Mann contains all the n-medium interactions followed by 
annihilation including antineutron rescattering in the initial state, the antineutron propagator 
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Gb is bare. The index "b" emphasizes this fact. Once the antineutron annihilation amphtude 
is defined by (8), the expression for the process amphtude (9) rigorously follows from (7). It 
is significant that the observable values (antineutron annihilation width F for example) are 
expressed through Mann- 

For solving the problem the field-theoretical approach with finite time interval is used [20] . 
It is infrared free. For the probability of the process (1) we get [11,16] 

Wb{t) ^ Wf{t) = eH\ Tt > 1, (10) 

where W/ is the free-space nn transition probability. Equation (10) does not contain density- 
dependence (F— dependence) since it corresponds to the limiting case Ft » 1 which is realized 
for the nn transition in nuclei. 

The lower limit on the free-space nn oscillation time is found to be 

r^in = 10'' yr. (11) 
This value is interpreted as the estimation from above. 



3.2 Model 2 with dressed propagator 

In the case considered above the amplitude Mann involves all the n-medium interactions followed 
by annihilation including the antineutron rescattering in the initial state. In principle, the part 
of this interaction can be included in the antineutron propagator. Then the antineutron self- 
energy E is generated and 

^'^^'^= e,-pi/2m-C/„-E = 
The process amplitude M^ is non-singular: 

Ma = eGaM'^^^. (13) 

The parameter S and " amplitude" of antineutron annihilation M^„„ are uncertain. The annihi- 
lation width F is expressed through Mann and not M^^„. Compared to Mann, ^ann truncated 
because the portion of the Hamiltonian "H is included in G^,- 

For the estimation we put M^^^ = Mann- This is an uncontrollable approximation. The 
process probability is found to be 

Wa{t) ^ Tat, 

^ ^r. (14) 
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The result differs from (10) fundamentally. Let r^j^ be the lower limit on the free-space nn 
oscillation time obtained by means of Eq. (14); Tnn is the oscillation time of neutron bound in 
a nucleus. Using condition Wd{Tnn) < 1, one obtains the relationships between t^^^ and Tnn: 

r^in^^^f^T- (15) 

We use the experimental bound on the neutron hfetime in oxygen Tnn > 1-77 • 10^^ yr obtained 
by Super-Kamiokande collaboration [8]. For estimation we take F = 100 MeV and E fa ReC4 — 
Un ~ 10 MeV {Un and Un are the potentials of n and n, respectively). Then 

= 1.2 • 10^ s, (16) 

which exceeds the lower limit given by the Grenoble reactor experiment [21] by a factor of 14 
and the restriction given by potential model [8] by a factor of five. Limit (16) is comparable to 
the restriction T^in(Fe) = (8 ^ 11) • 10^ s given in [15], however this coincidence is accidental. 
The model 1 should be rejected. If E — > 0, Wd rises quadratically. So r^j^ can be considered 
as the estimation from below. 

In the models with bare and dressed propagators the interaction Hamiltonians T-Lj and un- 
perturbed Hamiltonians are the same. The sole physical distinction between above-mentioned 
models is the zero antineutron self-energy in the model with bare propagator; or, similarly, 
the definition of antineutron annihilation amplitude {Mann or M^„„). However, it leads to the 
fundamentally different results. The problem is extremely sensitive to the value of antineutron 
self-energy E (see (9), (10) and (13), (14)) as well as the description of initial neutron state 
(propagator in the model 1 or wave function in the model 2) and the value of momentum 
transferred in the nn transition vertex [18]. This is because the amplitude (9) is in the peculiar 
point. 

On the one hand, the value E = 0, i. e. the model with bare propagator, seems quite realistic 
for pure physical reasons [11,18], and on the other hand the result is extremely sensitive to the 
details of the model. Besides, the calculation corresponding to the model 2 with bare propagator 
contains too many new elements and we view the result (10) with certain caution. Due to this 
we adduce the wide range of permissible values of T min : 

"^min '^min > Tmin- (17) 

This is a problem of great nicety. Further investigations are desirable. 

For the neutrons in bound state (see Fig. 2b) the results are the same as for nuclear matter 
[19]. 
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4 Summary and conclusion 



Since the operator (2) acts on the neutron, in the model 1 the vertex of virtual decay A — )■ 
n + {A — 1) is introduced because one should separate out the neutron state. This scheme 
is artificial because in the interaction Hamiltonian there is no term which induces the virtual 
decay A ^ n + {A — 1). The neutron of the nucleus is in the bound state and so it cannot be 
described by the propagator. In more exact terms, this is a crude model which is inapplicable for 
the problem under study. (This point is discussed more comprehensively in [22].) Alternative 
method is given by the model 2 which does not contain the above-mentioned vertex. The \in >- 
states are the eigenfunctions of unperturbed Hamiltonian what has no need of a commentary. 

It should be emphasized that oscillations in the vacuum and gas, in particular oscillations of 
ultra cold neutrons [6] are considered as well. These processes are described in the framework 
of the model shown in Fig. 2a [11,18]. They cannot be reproduced by means of model 1 in 
principle. For the reasons given above the model 1 should be rejected. If E — > 0, Wd{t) diverges 
quadratically. This circumstance should be clarified; otherwise the model 2 can be rejected. 
The calculation in the framework of the model with bare propagator gives the finite result, 
which justifies our approach from a conceptual point of view and consideration of the model 
with bare propagator at least as the limiting case. In reality this model seems quite realistic in 
itself [11,18,19]. 

Finally, taking into account the result sensitivity to the description of neutron state, we 
argue that the model 1 is inapplicable to the problem under study. The same is true for the 
potential model (see section 1). The model 2 and in the first place the value of antineutron self- 
energy invites further investigation. At present one can adduce only wide range of permissible 
values of t^iq. If E changes in the fimits 10 MeV > E > MeV, the lower limit on the 
free-space nn oscillation time is in the range 

10^^ yr > r^in > 1.2 • 10^ s. (18) 

Recall that for the free-space ab oscillations the ab transition probability is extremely sen- 
sitive to the difference of masses — rrib as well. The fact that process amplitude is in the 
peculiar point is the basic reason why the problem is extremely sensitive to E and the range 
of permissible values of Tmin is very wide. The values Tmin = 1.2 • 10^ s and Tmin = 10^^ yr are 
interpreted as the estimations from below (conservative limit) and from above, respectively. 
The estimation from below exceeds the restriction given by the Grenoble reactor experiment 
[21] by a factor of 14 and the lower fimit given by potential model [8] by a factor of 5. The 
range of uncertainty of Tmin is too wide. Further theoretical and experimental investigations 
are desirable. 
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